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ABSTRACT 

Fast radio bursts (FRBs) are extremely powerful sources of radio waves observed at cosmo- 
logical distances. We use a sophisticated model of FRB observations — presented in detail in 
a companion paper — to fit FRB population parameters using large samples of FRBs detected 
by ASKAP and Parkes, including seven sources with confirmed host galaxies. Our fitted pa- 
rameters demonstrate that the FRB population evolves with redshift in a manner consistent 
with, or faster than, the star-formation rate (SFR), ruling out a non-evolving population at 
99.9% C.L. Our estimated maximum FRB energy is logj Emax [erg] = 41.847097? (68% C.L.) 


-0.18 
assuming a 1 GHz emission bandwidth, with slope of the cumulative luminosity distribution 
y= -11589 15. We find a log-mean host DM contribution of 14516 pcem ? on top of a 


typical local (ISM and halo) contribution of ~ 80 pc cm^?, which is higher than most literature 
values. These results are consistent with the model of FRBs arising as the high-energy limit 
of magnetar bursts, but allow for FRB progenitors that evolve faster than the SFR. 
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1 INTRODUCTION 


Fast radio bursts (FRBs) are extragalactic transient radio sources of 
millisecond duration (Lorimer et al. 2007; Thornton et al. 2013). 
Some repeat, while most have not been observed to do so (Spitler 
et al. 2016; CHIME/FRB Collaboration et al. 2019b; Fonseca et al. 
2020a,b; Kumar et al. 2019; Shannon et al. 2018; James et al. 
2020), and the question of whether or not there are one, two, or 
more FRB populations remains open. The recent observation of a 
Galactic magnetar flare with FRB-like properties strongly suggests 
such objects as an FRB progenitor class (Mereghetti et al. 2020; 
The Chime/Frb Collaboration et al. 2020; Bochenek et al. 2020), 
though many more have been proposed (Platts et al. 2019). Yet, this 
flare was three orders of magnitude less powerful than the weakest 
FRBs, which in turn are orders of magnitude weaker than the most 
powerful FRBs (Shannon et al. 2018). FRBs may therefore have an 
unrelated origin. 

If the FRB population does originate from young magnetars, 
they would be expected to be closely correlated with star-forming 
activities, as observed for two rapid repeaters (Tendulkar et al. 2017; 
Marcote et al. 2020). However, the single largest sample of localised 
FRBs comes from the Australian Square Kilometre Array Pathfinder 
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(ASKAP; Bannister et al. 2019; Prochaska et al. 2019; Bhandari 
et al. 2020b). The host galaxies of these FRBs — which due to 
ASKAP's large field of view (FOV) and higher detection threshold 
tend to be the intrinsically most powerful bursts — do not show 
evidence for unusual star-forming activity (Bhandari et al. 2020a; 
Heintz et al. 2020). This allows for the possibility of much of this 
population to arise from other sources, e.g. compact binary mergers 
(see Caleb et al. 2018, and references contained therein). 


A useful method to distinguish between these models comes 
from the evolution of the FRB population on cosmological 
timescales. If FRBs originate from young magnetars, they will 
closely follow the star-formation rate (SFR) (Metzger et al. 2017), 
and peak in the redshift range 1—3. A binary merger scenario how- 
ever would likely lag the SFR, and possibly result in an FRB rate that 
is increasing with cosmological time (Cao et al. 2018). As yet, FRB 
population analysis has not been able to distinguish between these 
scenarios (Arcus et al. 2021; Luo et al. 2020). Other methods yield 
mixed results: Hashimoto et al. (2020) find evidence against the 
redshift evolution of once-off FRBs, and some evidence for redshift 
evolution of the event rate for repeating FRBs; while Locatelli et al. 
(2019) find evidence for an evolving FRB population for once-off 
FRBs. However, neither work follows a comprehensive approach 
advocated by Connor (2019), by modelling observational biases, 
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and allowing for the confounding effects of the FRB luminosity 
function. 

The FRB luminosity function is interesting in and of itself. 
Comparisons of the luminosity function of individual repeaters (e.g. 
Law et al. 2017) to the population as a whole tests the credibility 
of the one-population model, while evidence for a minimum burst 
energy above that produced by Galactic magnetars would require a 
separate progenitor class, or at least a separate emission mechanism. 
Models requiring rare events to explain FRBs can be challenged by 
measurements of the absolute volumetric rate (Ravi 2019). Esti- 
mates of the maximum FRB energy not only challenges theoretical 
models and pushes up against theoretical limits (Lu & Kumar 2018), 
but affects the ability to use FRBs as cosmological probes. Estimates 
of the host contribution to dispersion measure (DM) inform us of 
the environment surrounding FRB progenitors. Consequently, sev- 
eral groups have begun modelling the FRB population in an attempt 
to derive these parameters, although the results and methods have 
been inconsistent (Caleb et al. 2016; Luo et al. 2018; Lu & Piro 
2019; Luo et al. 2020; Arcus et al. 2021; Gardenier et al. 2020). 

In a companion paper (James et al. 2021), we present our 
method to model the FRB population. It uses the the methodology 
advocated by Connor (2019), and first implemented by Luo et al. 
(2020), while making several significant advances in accuracy and 
precision, and taking advantage of recent FRB localisations, and fit- 
ting for the measured signal-to-noise ratio. This models all known 
observational biases in detail, allowing us to make accurate and 
precise estimates of FRB population parameters, and model its cos- 
mological source evolution. Here, we present maximum-likelihood 
estimates of FRB population parameters using FRBs observed by 
the Australian Square Kilometre Array Pathfinder (ASKAP) and 
Parkes, and discuss the implications for the FRB population. 


2 REVIEW OF THE MODEL 


In modelling FRB observations, it is critically important to account 
for a range of observational biases. Our full treatment is contained 
in a (much lengthier) companion paper, James et al. (2021). To 
briefly summarise, we account for telescope beamshape, and re- 
duced observational sensitivity to high-DM, high-width FRBs, as 
recommended by Connor (2019); and fluctuations in cosmological 
dispersion measure according to best-fit cosmological parameters, 
local contributions from the Milky Way’s interstellar medium (ISM) 
and halo, and a log-normal distribution P(DMy, ost) of the host DM 
contribution, as per Macquart et al. (2020). This latter contribution, 
defined by 


2 
= (log DMp ost —Hhost)~ 


l ict (1) 


e 
y 
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is fit using the parameters Hhost and Chost- The effective host DM, 
DMhost; corrects the host DM for redshift: DMpogt = DM;, ost/ (142): 

Our model for the FRB population uses a power-law with 
cumulative slope y and maximum energy Emax, such that the prob- 
ability of observing an FRB above an energy threshold Ep is given 


by 
En Y Emax Y 
E E ) Emin = Enin 2 
p(E» Eh) = - l Enx y : (2) 


min 


P(DM host) = 


The minimum FRB energy is not well-constrained by current obser- 
vations, and is set to a very low value of 1030 erg. We scale the FRB 


energy E according to E ~ v®; for data taken exclusively at L-band 
(^ 1.4 GHZ), the model is almost degenerate to œ (a conclusion also 
reached by Lu & Piro 2019; Arcus et al. 2021), and so we use a 
symmetric Gaussian prior of a = —1.5 x 0.3 (Macquart et al. 2019). 

We model the evolution of the FRB population ®(z) (bursts 
per proper time per comoving volume) by smoothly scaling the SFR 
with the parameter n, 


= ®© (SR |" 
Pe), = Seal (63) 
We take SFR (z) from Madau & Dickinson (2014), 
2.7 
SFR(z) = li EA (4) 
os) 


Thus our full model treats Emax, y. Œ, n, Mhost aNd Chost as free 
parameters. 

There is some ambiguity in the interpretation of œ, which can 
instead be interpreted as a frequency-dependent rate. This is moti- 
vated by the many FRBs with a limited band occupancy, as origi- 
nally noted by Law et al. (2017) for FRB 121102. The interpretation 
of a; has slight effects on the modelling: in the rate interpretation, 
the FRB rate at high z is modified directly through a further factor 
@(z) ~ (1 +z)”, while in the spectral-index interpretation, this 
occurs indirectly, through the k-correction affecting the calculation 
of threshold energy Ej, from a fluence threshold Fin, and hence the 
rate via Eq. 2. In the absence of an obvious correct treatment, we 
by default present results for the spectral index interpretation, and 
show results for the rate interpretation in Appendix A. 

We use a sample of 24 non-localised, and seven localised, 
FRBs detected by ASKAP, and 20 FRBs detected by the Parkes 
multibeam system. These have been selected due to them occurring 
at high Galactic latitudes where the reduced sensitivity due to high 
Galactic DM is unimportant. The full telescope beamshape of each 
of these instruments is modelled in detail in our companion paper, 
based off the methods of James et al. (20192), while the reduction 
in sensitivity to high DMs and widths is modelled using the time- 
and frequency resolutions of the instruments according to Cordes 
& McLaughlin (2003). 


3 RESULTS 


Our single-parameter constraints are given in Figure 1, showing 
results both with and without a prior on a. Best-fit values and 
confidence limits, calculated using Wilks’ theorem (Wilks 1962), 
are given in Table Al. Two-parameter plots are given in Figure Al. 
We discuss the implications for each parameter individually below. 


3.1 Maximum burst energy Emax 


We find the maximum FRB to be energy is logo Emax (erg)= 
41 Bm (68% C.L.). According to our method, Emax is nor- 
malised to a bandwidth of 1 GHz at the mean frequency of the data 
used as inputs (about 1350 MHZ), and applies to all burst widths. 
A strict lower limit on Emax is set by the intrinsically brightest 
localised FRB, 190711, which — using a fluence of 34 Jy ms (Mac- 
quart et al. 2020), 1 GHz bandwidth, and œ = —1.5 — had an energy 
of E199711 = 10415 erg. Emax = 1041-6 erg is however consistent 
with observations at all quoted levels of confidence. 

The preferred value of Emax is most strongly correlated with a, 
which effectively attenuates FRBs as a function of redshift. Upper 
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Figure 1. Maximum likelihoods as a function of each considered variable (Emax, Œ, Y, n, Most; Chost) When marginalised over the other five, both with 
(orange, lower) and without (blue, upper) a prior on the spectral index a. Calculation results are given by points, with lines drawn using cubic spline smoothing. 
Vertical lines are single-parameter intervals at the labelled degree of confidence calculated using Wilks’ theorem with one degree of freedom. In the case of 


logy Emax, 90% and 95% lower limits are at 41.4. 


limits on Emax are also strongly correlated with y, since a large 
negative value of this parameter makes it unlikely to observe FRBs 
near Emax- 

Our value of Emax lies in the middle of the values found by 
other authors. From Figure A1, fixing n = 0 as per Luo et al. (2020) 
would lead to a lower value of Emax, and greater consistency with 
that work. The higher values of Emax found by Lu & Piro (2019), 
and used by Arcus et al. (2021), arise in models that assume a 1-1 
DM-z relation, which will tend to over-estimate Emax when an FRB 
with a significant excess DM — either due to its host or intervening 
matter — is detected. 

A key implication of Emax is the distance out to which an 
FRB is observable by a given telescope. For œ = —1.5, our value of 
logjo Emax leads to a maximum observable redshift of z = 4*2 as 
for an instrument with 1 Jy ms threshold. 


3.2 Intrinsic luminosity index y 


Our best-fit power-law index for the FRB population is y = 
-]. i er (68% C.L.). As discussed by Macquart & Ekers (2018), 
this parameter primarily governs the degree to which FRBs are 
viewed from the near or far Universe, with steep values of y (i.e. be- 
low -1.5) leading to observations being dominated by nearby events 
and the event rate being governed by Emin. Our result is definitely 
above this value, which is in agreement with all other calculations. It 
is however somewhat steeper than the values found by other authors. 

Why? Luo et al. (2020) assume no cosmological source evolu- 
tion, which this parameter is strongly correlated with. An increase of 
high-redshift FRBs can be due to either a lower y, leading to more 
bursts visible near Emax in the larger volume of the distant Uni- 
verse; or due to an evolving population, as determined by n. This 
anti-correlation is clearly visible in Figure Al. Both Arcus et al. 
(2021) and Lu & Piro (2019) allow source evolution, but assume a 
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1-1 DM-z relation, implying a large distance for the highest-DM 
FRBs. In order to fit such bursts without over-predicting a large 
number of lower-energy bursts requires a flat luminosity function. 

In the case that all FRBs repeat, with each FRB having the 
same Emax and y but a distribution of intrinsic rates, the intrinsic 
luminosity function for the entire population will match that of each 
FRB. This index has been well-measured for FRB 121102, with 
data giving a range y121102 ~ —0.9 + 0.2 (Law et al. 2017; Gajjar 
et al. 2018; James 2019). This is consistent with our value for the 
population. However, should Emax vary over FRBs (which is quite 
likely), then the value of y for the population might be steeper. 
FRB observations by CHIME, which have detected several FRBs 
with many repeat bursts (CHIME/FRB Collaboration et al. 2019b; 
Fonseca et al. 2020a,b), should be able to answer the question 
definitively. 


3.3 Redshift evolution n 


Our best-fit value of the redshift evolution scaling parameter is 


n= De (68% C.L.). Under the ‘rate interpretation’ of a, we 


find n = 1.26*031. In both cases, n = 0 is excluded at better than 
99.9%, which holds when no prior on a is considered. 

Our detection of evolution in the FRB population supports 
conclusions based on FRB localisations, which locate most FRBs 
within normal host galaxies (Heintz et al. 2020); evidence associ- 
ating FRBs with magnetars, such as the recent Galactic magnetar 
outburst (The Chime/Frb Collaboration et al. 2020; Bochenek et al. 
2020); and observations of the host environment of FRB 121102 
(Michilli et al. 2018), as well as predictions from several classes of 
progenitor models (Platts et al. 2019). 

This does not mean that we have confirmed that FRBs exhibit 
cosmological evolution identical to the star-formation rate however. 
A more-general model of source evolution, as used by Lu & Piro 
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Figure 2. Maximum likelihood fits for source evolution parameter n, for 
three different cases: interpreting a as a spectral index, and source evolution 
n scaling the star-formation rate as per Eq. 4 (blue); interpreting @ as a 
frequency dependent rate, again using Eq. 4 for n-scaling (orange); and with 
a as a frequency-dependent rate, but n scaling source evolution simply as 
(1 + z)?-”" (green). The vertical lines show 90 C.L. intervals calculated 
using Wilks' theorem. 


(2019), simply assumes a (1 + z)" dependence, i.e. it removes 
the denominator and normalising constant in Eq. 4. Near z = 0, 
n' = 2.7n — however, the models will diverge above z = 1. A 
true detection of scaling with the star-formation rate would require 
observations to be consistent with a downturn relative to the ( 12) 
model at and beyond the peak of star-forming activity. 

Figure 2 plots the likelihood for both interpretations of a, and 
the (1 + z)?-7" model. While the spectral index interpretation of 
a clearly gives a better fit, the difference in maximum likelihoods 
between the two source evolution models under the rate assumption 
is negligible, with the preferred value of n being slightly higher 
under SFR-scaling to compensate for the denominator in Eq. 4. 
We therefore cannot claim evidence for a down-turn in the source 
evolution function due to the peak in the SFR, only that the FRB 
population is evolving with cosmological time, with the rate per 
comoving volume greater at higher z. 

Our result is still a significant improvement on prior works. 
Previous calculations have either had to assume a value for FRB 
source evolution (of n = 0 or 1), due to complete degeneracy with 
a (Lu & Piro 2019), or otherwise could not distinguish between 
models (Caleb et al. 2016; Lu & Piro 2019; Arcus et al. 2021). As 
previously noted, and explained in detail in our companion paper, 
this degeneracy also affects this work. However, the degeneracy is 
not complete — it is partially broken by the ASKAP/ICS sample 
of localised FRBs, and by the Parkes sample, which probes to suffi- 
cient z to be sensitive to the non-Cartesian nature of the Universe. 
Of similar works, only Caleb et al. (2016); Luo et al. (2020) model 
a telescope beamshape. We show in our companion paper that these 
authors’ assumption of a Gaussian (~Airy) beamshape for Parkes 
observations is sufficient, but doing so for ASKAP data — as con- 
sidered by Luo et al. (2020) — is inappropriate. The inability of 
Arcus et al. (2021) and Lu & Piro (2019) to exclude n = 0 may 
be due to their lack of beamshape modelling. Including beamshape 
reveals that a larger fraction of the sky is probed at lower sensitivity, 
thus increasing sensitivity to FRBs in the local Universe relative to 


that in the distant Universe. Without this effect, n must be artificially 
decreased to model the observed number of near-Universe bursts. 


3.4 Excess DM distribution 


Our model fits a log-normal distribution to DMhost, which nomi- 
nally covers the host galaxy and the immediate FRB environs. The fit 
will naturally include deviations from the NE2001 DM model of the 
Milky Way and the assumed halo DM of 50 pc cm ?. We find best- 
fit values of logi Hhost = 2107791 and logio Chost = 0:5120, 
with both parameters being relatively independent of the other 
four. Figure Al shows that high values of post and low values 
of Ghost are most strongly excluded. The only other authors to 
fit these parameters are Macquart et al. (2020), who use a sub- 
set of the data analysed in this work; our fitted value for the 
mean DM is greater than theirs, but not significantly. Partially, 
this is because Macquart et al. (2020) do not account for re- 
duced sensitivity to high-DM bursts. Our inclusion of this effect 
requires a greater intrinsic high-DM population to fit the same 
observations. Combined with local contributions from the Milky 
Way’s ISM and halo, we estimate a mean non-cosmological DM 
of DMysm + DMhalo + Hhost = 50 + 35 + 145 = 230pccm™? at 
z = 0. However, this still allows for low values of DM observed 
by ASKAP (Shannon et al. 2018) and CHIME (CHIME/FRB Col- 
laboration et al. 2019a), since both DMpog and the cosmological 
contribution can vary. This large value of mean non-cosmological 
contribution helps to explain the observation by Shannon et al. 
(2018) that the mean DM of the Parkes FRB sample is not as large 
relative to the ASKAP/ICS sample as would be expected from the 
relative telescope sensitivities alone. 


3.5 The prevalence of FRBs 


We estimate the best-fit absolute rate of FRBs above Emin; Po, 
by maximising the product of p, between the ASKAP/FE and 
Parkes/Mb samples with well-constrained Tops. We quote 39, de- 
fined as the estimated rate of bursts above 103° erg (above the max- 
imum allowed value of Emin James et al. 2021) per year at z = 0. 
In the case of our best-fit model, we find 39 = HEY - 10^ bursts 
Gpc ? yr-! (90% C.L.). 

This value is broadly consistent with that estimated by other 
authors (Ravi et al. 2019; Lu & Piro 2019; Luo et al. 2020), and 
supports the conclusion that the majority of FRBs must either be 
repeaters, or cannot be due to known populations of once-off events. 

Interestingly, the best-fit parameter set under-predicts the num- 
ber of FRBs observed by ASKAP/FE (12.9 vs. 20 in 1274.6 days), 
and over-predicts the number found by Parkes (17.0 vs 12 in 164.4 
days). There are several possible causes of this discrepancy. One 
possibility is a minimum FRB energy — or at least a flattening of 
the distribution at low energies — which would reduce the number 
of bursts seen by the more-sensitive Parkes telescope. Another is 
the low number of FRBs detected by Parkes with SNR below 16, 
as noted by James et al. (2019b), which could be an indicator of 
a reduced detection efficiency to low-fluence bursts. Both are in- 
vestigated and deemed unlikely in our companion paper. A third 
option is that the observation times reported here are raw observa- 
tion times, and do not account for lost effective observation time due 
to e.g. radio-frequency interference, which is likely more prevalent 
at Parkes than ASKAP. Finally, this could be the result of simple 
statistical fluctuations — if the true rates are 12.9 and 17.0 FRBs 
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respectively, then the product of the ASKAP and Parkes likelihoods 
will be this or less unlikely 7.7% of the time. 


4 CONCLUSION 


We have used a precise and accurate method of modelling the results 
of FRB surveys to fit the measured DM, z, and signal-to-noise 
ratios of FRBs detected by ASKAP and Parkes. We have carefully 
selected our data to ensure it is not biased due to under-reporting of 
observation time, or due to large local DM contributions reducing 
sensitivity. Crucially, we have included a sample of localized FRBs 
from ASKAP for which the redshift of the host galaxies is measured. 

These modelled observations are tested against a six-parameter 
model of the FRB population. Using a maximum-likelihood ap- 
proach, we have derived the tightest constraints on FRB popula- 
tion parameters to date. Our value of the maximum FRB energy 
of 41.845009 erg (68% C.L.) is mid-way between previous esti- 
mates.The intrinsic slope of the cumulative luminosity distribution, 
y, is found to be —1.16*0-1 (68% C.L.), consistent with, but slightly 
steeper than, the slope found for FRB 121102. Importantly, we find 
that the FRB population evolves with redshift, scaling with the 
star-formation rate (SFR) to the power of Lg or 1207 3. 
depending on the interpretation of FRB spectral properties. While 
we cannot distinguish between SFR-scaling and a model where the 
FRB population increases as a simple power of ( 14-2)?" , in all sce- 
narios we exclude a non-evolving population at better than 99.9% 
CL: 

Our best-fit log-mean host contribution to DM of 145 pc cm™ 
is also somewhat higher than the standard value of 100 pc cm" ?. 

Such large excess dispersion measures, and a population evo- 
lution consistent with star formation, strongly aligns with the hy- 
pothesis of FRBs originating from young magnetars. We caution 
that these results apply to the total FRB population (which may or 
may not consist of multiple sub-populations), and only to that part 
of the population to which the ASKAP and Parkes observations are 
sensitive. 
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APPENDIX A: FURTHER DATA 


Figure A1 shows the correlation plots between all two-parameter 
combinations, excluding o. 

In the main body of this work, we present parameter limits at 
90% C.L. obtained when using « as a prior. Table Al presents pa- 
rameter limits are different confidence levels, both with and without 
a prior on a. Furthermore, in Section 2, we note that the spectral 
index a can also be interpreted as a frequency-dependent rate. Our 
standard set of results use the spectral index interpretation. For com- 
pleteness, in Table A2, we give the single-parameter limits under 
the rate interpretation. 


This paper has been typeset from a TEX/IATEX file prepared by the author. 


MNRAS 000, 1—6 (2020) 


FRB source evolution 7 


-2213 
E 
p 
-222 8 
9 
-223 5 
5 
+ 
—224 
-220 
^x 
—221 £ -221 _ 
E S = 
-2228 -222 6 E 
E E à 
-223 2 -223 Š 8 
i E E 
A 
-224 Ẹ 
-220 -220 
-221 _ -221 
3 c 
-222 8 -222 8 
i Š 
> $ 
-223 8 á 
Ë -223 f 
-224 E 
-220 -220 
-221 -2214 
z E 
-222 3 -222 i 
5 2 
È 8 
-223 f -223% 
+ E 
- 
-224 -224 


Figure A1. Two-parameter maximum likelihood results, showing 68%, 90%, and 95% confidence intervals, calculated using Wilks’ theorem and a x3 
distribution. 
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No prior Prior on « 
Parameter BestFit 68% C.L. 9090 C.L.  959c C.L.  BestFit 68%CL. 9096 C.L. 95% C.L. 
Gie a -- mr 
a -2.50 N/A N/A N/A -1.55 N ecd er 
y 445 — 7$ NAT a 46 — "wb "s ME 
ý 2.22 "S “0.96 Mis 17 0 "Ws “0.66 MES 
Hhos 214 —— T" “Oa "us 2106 — "$5 NE NET 
Fhost 0534 — T ^is "Qm 09 — "p Myr d 
Table A1. Confidence limits on single parameters, both with (left) and without (right) a prior on a. 
No prior Prior on « 
Parameter BestFit 68% C.L. 9090 C.L.  959c C.L.  BestFit 68% C.L. 9096 C.L. 95% C.L. 
hg E« 45 22 — 32 — 38. «x ey 33 — 33 
a -2.50 N/A N/A N/A -1.50 oU 1022 be 
y 449 — "ui "as uno 5 X 55 — X 
AEG uu 26 Wl a M 
B aw uu WW WW M 0 M 
g dE uu NM M M M M 


Table A2. Equivalent of Table A1, calculated assuming the rate interpretation of œ (see Section 2). 
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